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1 .   INTRODUCTION 

The  purpose  of  this  paper  is  to  investigate  a  general  plan  for 
the  local  recognition  of  visual  patterns.  This  local  recognition  is  the 
basis  of  a  higher  level  recognition,  global  recognition  [k] . 

The  following  subjects  are  investigated  in  our  local  recognition 
level: 

1.  A  General  Labeling  Routine, 

2.  Cleaning  of  Random  Noise, 

3.  Shaping,  including  Burst  Noise  Correction,  and 
k.  Preliminary  Extraction  of  Feature  Objects  (the 

basis  of  higher  level  recognition). 

The  local  filtering  theory  follows  lines  developed  by  Yamada  [l]. 
Programs  to  test  the  filtering  theory  were  written  in  PAX  [2],  a  parallel 
processing  language  written  for  the  IBM  709^-1^-01  System  at  the  University 
of  Illinois.   Portions  have  been  written  in  SCATRE  [3],  the  assembly  language 
for  that  system.  An  additional  routine  not  described  in  the  original  PAX 
report  has  been  used.   This  routine  is  described  in  Appendix  I. 


[l]   Yamada,  S.,  "Local  Recognition  of  Pictorial  Images,"  Department  of 

Computer  Science,  University  of  Illinois,  Urbana,  Illinois.   Report  in 
preparation. 

[2]   Stein,  J.  H.,  "User's  Manual  for  PAX,"  Department  of  Computer  Science, 
University  of  Illinois,  Urbana,  Illinois.   Report  No.  IU7  (July,  1963). 

[3]  Wolf,  W.,  "SCATRE  for  the  IBM  709U.   Library  Subroutine:  L1-U0I-SCRE-19-BX, " 
Department  of  Computer  Science,  University  of  Illinois,  Urbana,  Illinois. 

[h]      Yamada,  S.,  "Global  Recognition  of  Pictorial  Images,"  Department  of 

Computer  Science,  University  of  Illinois,  Urbana,  Illinois.   Report  in 
preparation. 
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2 .   LABELING 

2.1  General 

A  subroutine  which  is  common  to  both  programs  in  Parts  3  and  ^  is 
known  as  LABEL .   As  the  name  implies,  the  routine  labels  the  points  of  a 
picture  according  to  the  approximate  direction  of  the  line  through  the  point 
The  labeling  convention  is  as  follows: 


Actually",  the  program  assigns  at  least  two  directions  to  each  point 
labeled.   If  all  lines  were  straight  lines,  only  one  label  would  be  necessary 
(i.e.,  a  line  running  east  and  west  would  need  only  one  label  (l)  rather  than 
two  labels  (l  and  5))«   However,  a  point  on  a  curved  portion  of  a  line  or  on  a 
corner  of  a  line  may  require  two  directions  which  are  not  diametrically  opposed 
(i.e.,  the  labeling  directions  may  be  1  and  k) . 


2.2  Use 

The  subroutine  is  entered  at  the  symbol  LABEL  using  PAX  index 
register  8  as  a  link.   The  picture  to  be  labeled  is  in  a  plane  called  INPUT. 
Upon  return  to  the  calling  program,  INFUT  is  unchanged,  all  labeled  points 
are  in  a  plane  called  NEWVAL  (new  values),  and  all  unlabeled  points  are  in  a 
plane  called  NOCON  (no  connections),  and  the  labels  are  in  eight  planes 
beginning  with  IODIR  (input-output  directions).   Plane  IODIR-1+k  contains  all 
points  labeled  in  the  k  direction.   (Note:   NOCON  =  NEWVAL)   Two  other  planes, 
TEMPI  and  TEMP2,  are  used  as  temporary  scratch  planes  by  the  subroutine . 
PAX  index  registers  1  through  h   are  altered. 
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2,3  Process 

The  initial  labeling  is  done  in  a  somewhat  unconventional  manner, 
Rather  than  look  at  immediate  neighbors,  the  program  examines  points  two 
units  away  along  the  primary  directions,   A  point  is  labeled  in  directions  d 
and  d+k   if  the  point  itself  and  the  second  neighbors  along  the  d  and  d+h 
directions  are  all  black  in  INPUT,   The  labeled  points  are  added  to  NEWAL, 
After  all  directions  have  been  examined,  NOCON  is  set  to  the  complement  of 
HEWVAL. 

After  this  initial  labeling,  the  IODIR  planes  should  contain  enough 
information  about  the  straight  portions  of  the  input  picture  that  the  corners 
and  gaps  can  easily  be  filled  in  without  referring  to  INPUT. 

The  filling  process  is  accomplished  by  examing  immediate  neighbors, 
A  point  is  labeled  in  the  d  direction  if  its  d  neighbor  and  either  its  d+3, 
d+U,  or  d+5  neighbor  has  been  labeled  in  the  initial  labeling.   This  allows 
points  to  be  labeled  along  straight  portions,  or  around  curves  or  corners „ 

Points  which  were  labeled  in  the  initial  labeling  are  not  allowed 
to  pick  up  additional  labels  in  the  filling,,   This  is  done  by  filling  only 
points  in  NOCON „   Note  that  filled  points  are  not  subtracted  from  NOCON  until 
all  filling  is  completed. 


2 ,h     Results 

This  labeling  algorithm  has  proven  to  be  useful  in  two  cases.   The 
first  is  where  the  lines  of  a  picture  have  a  width  of  from  one  to  three  points 
and  are  very  long  compared  with  their  width.   When  wider  lines  are  used,  the 
interior  points  are  multiply  labled,  but  the  border  points  of  the  line  are 
labeled  according  to  the  direction  of  the  border.  Therefore,  this  routine 
performs  point  labeling,  rather  than  line  labeling. 

The  second  case  is  where  the  lines  are  embedded  in  random  noise 
(see  Section  3), 

The  most  prominent  fault  of  this  routine  is  that  the  lines  are 
''chopped,"   Since  the  labeling  depends  on  the  continuity  of  the  line  through 
the  point,  the  end  points  of  the  line  are  not  labeled.   In  the  case  where 
noise  is  not  prevalent,  it  is  a  simple  matter  to  extend  the  labels  to  the  ends 


1,   See  program  listing  in  Appendix  I  and  flow  chart  in  Figure  1, 


of  the  INPUT  lines.   This  is  done  in  the  labeling  part  of  the  reconnection 
routine  (Part  k) .      However,  where  noise  is  prevalent,  the  extension  becomes 
extremely  difficult,  and  therefore  was  not  used  in  the  noise  cleaning 
routine  (Part  3). 


Enter: 
LABEL 


I 


IODIR(d)  =  Black  points  in  INPUT  with  second 
neighbors  in  d  and  d+^  directions  marked 


i. 


NEWVAL  =  IlODIR(d) 

d 

NO CON  =  NEWVAL 


t 


NEWPTS  =  NOCON  •  IODIR(d)  •  [IODIR(-d)  +  IODIR(-d+l)  +  IODIR(-d-l)] 
IODIR(d)  =  IODIR(d)  +  NEWPTS 


NOTE:   "-d"  = 


d+U  for  d<l+ 
d-h   for  d>i+ 


NEWVAL  =  IlODIR(d) 

d 

NOCON  =  NEWVAL 

} 

r 

RETURN  to  Main  Program 

Figure  1.   LABEL  Subroutine  Flow  Chart 
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3  .   NOISE  CLEANING 

3.1  General 

This  noise  cleaning  routine  is  a  simulation  of  a  general  algorithm 
worked  out  by  Yamada  [l]  for  special  application  to  large  noise  cases.   The 
large  noise  case  was  chosen  since  it  is  a  more  general  problem.   If  no  a  priori 
information  is  known  about  the  extent  of  noise  in  a  picture,  the  most   power- 
ful cleaning  routine  must  be  used.   On  the  other  hand,  if  it  is  known  that  the 
noise  is  small,  a  less  powerful  but  more  efficient  routine  may  be  used. 

This  cleaning  algorithm  cannot  be  programmed  precisely  because  it  is 
neither  definite  enough  nor  simple  enough  to  be  programmed  with  any  degree  of 
efficiency.   Also,  the  algorithm  is  based  on  a  well-defined  noise  input.   As 
will  be  seen  in   paragraph  3»2,   the  noise  used  is  not  well-defined. 

However,  the  original  filtering  theory  has  been  modified  and  realized  as 
a  programmable  and  fairly  efficient  routine..   The  label  routine  described  in 
a  previous  section  uses  aspects  of  the  general  theory.   This  section  describes  the 
extension  of  the  label  routine  into  a  noise  cleaning  routine, 


3 .. 2  Generation  of  Random  Noise 

A  pattern  was  chosen  as  a  basis  of  study,  and  noise  was  added  to 
this  pattern.   The  noise  was  to  be  of  a  random  nature,  but  of  uniform 
density  over  the  entire  plane. 

In  an  effort  to  achieve  this,  a  random  number  generator  was  used 
to  generate  random  coordinates.   The  random  generator  routine  was  entered 
twice,  and  the  two  numbers  obtained  were  used  as  x-  and  y-  coordinates  of  a 
point,.   That  point  was  then  made  "bad"  by  reversing  its  value.   If  the  point 
had  previously  been  reversed,  no  change  was  made, 

The  density  of  noise  was  controlled  by  the  number  of  points  which 
were  changed.   By  using  a  noise  level  scale  from  0  to  9,    the  entire  range  of 
possible  degrees  of  noise  was  covered,  A  level  of  n  indicates  that,  on  the 
average,  n  points  of  each  3x3  window  are  "bad,"   In  other  words,  in  an 

Q 

m  x  m  plane,   nm  /9  points  were  made  bad  in  order  to  achieve  a  noise  level 
of  n. 
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The  noise  obtained  by  this  generation  routine  is  usable,  but  far 
from  ideal.   There  are  two  reasons  for  this.  First,  the  random  number 
generator  did  not  generate  uniformly  from  0  to  71,  but  put  a  slight  emphasis 
on  values  from  0  to  5^.   Second,  it  was  found  that  sets  of  two  randomly 
generated  numbers  do  not  define  random  coordinates  in  a  two-dimensional  plane. 

3 .3  Noise  Cleaning  Process 

The  basis  of  the  program  is  the  LABEL  routine.   But  before  attempting 
to  label  a  noisy  picture,  an  initial  mapping  of  the  domain  to  the  point  is 
performed.   That  is,  only  points  (black  or  white)  with  greater  than  or  equal 
to  a  certain  number  of  black  neighbors  (including  itself)  are  input  as  black 
points  to  the  label  routine.   This  number,  T,  is  called  the  threshold  level. 

Because  of  the  nonuniformity  of  the  noise,  no  single  T-value  is 

optimum,   Using  T=U,  too  many  bad  points  are  retained_,  and  the  picture  is 

not  thoroughly  cleaned.   Using  T=5,  not  enough  points  are  retained  and  small 
portions  of  lines  are  erased. 

Therefore,  the  picture  to  be  cleaned  is  labeled  twice,  using  T=^  and 
T=5 o   The  portions  erased  by  using  T=5  are  picked  up  by  extending  the  results 
of  T=5  into  the  results  of  T=U  on  a  direction  continuation  basis.   The 
resulting  picture  is  then  reextended  into  T-h   until  no  new  points  are 
picked  up. 

The  extension  is  begun  with  the  results  of  the  labeling  using  T=U 
in  four  planes  starting  with  Ik,      T^(d)  contains  the  points  labeled  in  the 
d  and  d+k   directions.   Similarly,  T5(d)  contains  the  d  and  d+k   labeled  points 
using  T=5„   T5ALL  contains  all  points  labeled  using  T=5.   T5BAR  =  T5ALL,  and 
contains  all  points  not  labeled. 

The  end  points  in  T5  are  selected  by  direction  and  are  put  into  one 
of  eight  planes,  DENDS(d).   A  point  is  an  end  point  in  the  d  direction  if  it 
is  in  T?(d),  has  no  labels  in  the  side  or  corner  in  the  d  direction,  and  has 
any  label  in  the  opposite  direction  'see  Example  l) .   In  an  effort  to  retain 
continuity  at  corners,  all  points  in  TU(d+l)  and  TU(d-l)  are  added  to  DENDS(d) 
to  form  DENDSX(d). 


1.   See  program  listing  in  Appendix  n  and  flow  chart  in  Figure  2, 
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Example  1.   Examples  of  End  Point  Definitions 


Requirements  for  end  points  in  the  2  direction: 


D 

N 

N 

D 

X 

N 

L 

D 

D 

Requirements  for  end  points  in  the  3  direction: 


N 

N 

N 

D 

X 

D 

D 

L 

D 

In  these  examples,  X  is  the  point  in  question  and  must  be  labeled 
in  the  same  direction  as  the  end  point  is  to  be  labeled;  N  is  a  point  with 
No  label;  L  is  a  Labeled  point  (regardless  of  direction);  and  D  is  a  "Don't 
Care"  point  (does  not  influence  the  selection  of  end  points). 
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Read  ORIG  picture 


INPUT  =  Points  in  ORIG  with  >  5  black  neighbors 


[label  input"] 


T5(d)  =  IODIR(d)  +  IODIR(d+U) 
T5ALL  =  NEWVAL 
T5BAR  =  NOCON 


INPUT  =  Points  in  ORIG  with  >  k   black  neighbors 


LABEL   INPUT 


I 


TU(d)    -  IODIR(d)    +  IODIR(d+U) 


DENDS(d)  =  End  Points  of  lines  in  T5(d) 


I 


DENDSX(d)  =  DENDS(d)  +  TU(d+l)  +  TU(d-l)         |  T|?(d)  =  T5(d)  +  NEWPTS(d) 


1 


NEWPTS(d)  =  Points  found  by  chaining 
DENDSX(d)  into  Th(d) 


A 


A 


COMP  =  NE¥PTS(d)    -    T5(d) 


Null 


jfcj   Print  T5's  | 


Figure  2.   Noise  Cleaning  Flow  Chart 
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The  growing  is  accomplished  by  chaining  the  points  in  DENDSX(d)  into 
TU(d).   The  new  points  are  put  into  NEWPTS(d). 

After  all  directions  have  been  grown,  NEWPTS(d)  is  tested  against 
T5(d).   (NEWPTS(d)  may  contain  points  previously  in  T5(d).)   If  no  new  points 
have  been  picked  up,  the  process  is  complete  and  the  T5 ' s  contain  the  final 
picture.   If  new  points  are  found,  they  are  added  to  the  T5's  and  the  extension 
is  reiterated. 


3.U  Results 

The  program  was  tested  using  noise  levels  from  0  to  3.   The  results 
are  shown  in  Figure  3-   Comparison  of  input  and  output  pictures  show  that  all 
prominent  features  of  the  input  are  retained  in  the  output. 

The  routine  obviously  does  not  work  well  in  the  highest  noise  case, 
but  this  degree  of  noise  is  extremely  rare  in  actual  pictures  produced  by 
electrical  or  mechanical  means. 

The  symbols  which  mark  points  on  the  processed  pictures  of  Figure  3 
specify  the  labels  assigned  to  the  points  in  accordance  with  the  code  of  the 
table  below.   For  example,  a  point  on  a  processed  picture  which  is  marked  "3" 
was  at  the  juncture  of  a  "2-6"  (Northeast-Southwest)  line  segment  and  a  "3-7" 
(North-South)  segment. 
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Direction 

Code 

Symbol 
.  Di  it. pi  it. 

on 
Pi- 

ed 

nrp 

1-5             2-6 

3-7 

k-Q 

E 

X 

A 

X 

N 

X 

B 

X 

1 

X                       X 

2 

X 

X 

3 

X 

X 

k 

X 

X 

6 

X 

X 

+ 

X                       X 

X 

- 

X                       X 

X 

* 

X 

X 

X 

• 

X 

X 

X 

$ 

X                       X 

X 

X 

>  1  direction 


>  2  directions 


*    3  directi 


ons 


k   directions 
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xxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxx 

xxxxxxxx        xxxxxxxx 

xxxx  xxxxx 

xxxxx  xxxx 

xxxx  xxxx 

XXX  XXX 

XXX  XX 

XXX  XXX 

XXX  XXX 

XX  XXX 

XXX  X 

XX 
XX 
XX 
XX 
XX 
XXX 
XX 

XXX     xxxxxxx 

XX   xxxxxxxxxxxxxxx 

xxxxxxxxx     xxxxx 

XXXXXXXX  xxx 

xxxxx  xxx 

XXX  XX  xxx 

XXX  XX  xxx 

XXXXX  XX 

XXX  XX  XX 

XX  XX  XX 

XX  XX  XX 

XX  XX  XX 

XXX  XX  XX 

XXX  XX  XX 

XX  XX  XX 

XX  XX  XX 

XXX  XX  xxx 

XX  xxx  xxx 

XX  XX  xxx 

XX  xxx  XX 

XX  Xa  xxx 

XX  xxxx  xxx 

xxx  xxx       xxxx 

XX  xxxx    xxxx 

XX  xxxxxxxxxx 

XX  xxxxxx 

XX 
XXX 
XX 
XX 
XXX 
XX 
XX 

xxx 

XX 
XX 

xxx 

XX 
XX 
XX 

xxx 

XX 
XX 

>xx 

XX 
XX 
XX 

xxx 

XX 
XXX 
XX 

xxx 


ORIGIN/SI    PICTURE  ERROR  LEVEL  =   0 


) 

ex 

xxxx 

xxxxx 

XX 

XX 

XX 

XX 

xxx 

XX 

XX 

XXX 

XXX 

xxx 

XX 

XX 

XX 

XX 

XX 

XXX 

XX 

XX 

XX 

XXX 

XXXXXX 

XX 

XXXXXXXXX     XX 

XXX     XXXXXX  xxx 

XX 

xxxxxxxx 

XX 

XXXXXX 

XX 

XX 

XX 

XXX 

XX 

XX 

XX 

XX 

XX 

XX 

X 
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EEEEEEEfc 

EEEEEEEEEEEbtEEEE 
1EEEEEE  Et^EC 

111  ^4^BB 

1111  ^44B 

All  B4B 

AA  BBB 

AA  BBB 

AA  BBE 

AA  B 


AAA 

AA 
AA 
AA 
A 
3 
3N 
33N 
NN6 

3NN       EEE 
33«*EEEEEEEEEEEL 
A32EEEtE         E44 
A1+11E  BB6B 

11E23A  RB6 

A1E2NA  B66 

AA   NN  N 

AA    NN  NN 

A  A     NN  NN 

3  A       NN  NN 

33       NN  NN 

33       NN  NN 

33        NN  NN 

33        NN  N 

33         NN  N 

NN         NN  A3 

N  N  33 

N  N  33 

NN  66  333 

NN  66B  A33 

NN  B6B  AAA 

N  eEB  AAA 

N  BAB  HA 

N  BB4       11A 

NN  ^EEEEEEl 

NN  EE 

NN 
N 

NN 
N 
6 

NN 

NN 

6 

N 

NN 
N 
N 

NN 
N 
N 
N 


N 

N 

NN 

NN 
6 

N6 

N6 

N 


PROCESSED    PICTURE  ERROR  LEVEL  *  0 


Al 

EB 

331466 

33 

B6 

33A 

B66 

3N 

N6 

33 

6 

N3 

NN 

33 

N6 

33N 

66 

3NN 

NN 

33N6 

No 

3  +  444 

36N 

331EE44E 

3N6 

3A     ^E4^EE+266 

3A 

EEEEL*6 

3 

E4466 

3 

BB6 

N 

B66 

B6 

N 
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X  X  XX 

X  X  X  X        X  X  XXX    XXXX     XXX  X 

XX  XX    XXXXXXXX    XXX    X    xxxxxx 

X  XX  XXXXXXXX       X  XXXX    XX     X 

X  XX       XXXX  XXXXXX        X 

X  XX  XXXXX  X  XX        X       XXXX     X        X 

X  XX     XXX  XX  XX  XXXX 

XX  X  XXX  X  X     XXX        X 

X  XXX  XXX     XX  XXX  XX 

X  XXXXX  X  XXX 

X  X   XXX     X        X       XX  XXX 

X  X   X      XX      XX  X  X        XXX 

X  XX  XXX  XX  X 

X   XX  XX    XX  X      XX 

X  X    X   XX  X  XX  X 

X  XX  XX 

XX  XX  XXX  X  XXXXX  XX 

X     X  XX  XXX  X 

XX  X     X        XX  XX 

XX  XX      X  X  X  X  XX 

X  X       XX  XX  XX   X  X  X 

X  X        X   XX  XX  XXXXXXXXX  X  XX  X 

X  X    XXXX  X  XX       XXXXXX       XX  X 

X  XXXXXXXXX  XXXX   XX  XX 

XXXXX  X  X       XXX  X         XX       X 

X   XXX   XX  XX  XXXX  X 

X  XXX   XX  X  X         XXX  X   X 

X  XXX    XXX  XX  X 

X         XXX      XX         X       X  X       XX  XXX 

X  XX       XXX  X        XX  XX  X 

X     XXX       XX   X      X        X       XX         X   XX      X  X  XX       X 
XX         XX        XX  XX  X 

X  XX  XX  XX  X     XXX  XXX  X 

XXXX        XX  XXX  XX  X  X 

X     XX       X     X  XX  X        XX       X     X  XX  XX 

XX  X  XX  X        X     XX  XX  X 

X     XXX  XXX  X    XXX  X  X 

XX    X  XX  X    XXX  XX  X  X 

XXX  XX  X    XXX  X 

XX  XX         XXX   X  X    X        XX        X  X 

X   XX        X      XX  X  X       XX  XX 

XX  XXXXXX  XXX         XX  X 

X       XX     XX         XXX   X      XX  XX    X  X  X 

XX   X       X      XXXX      XX  X  X       XX       X  X 

X   XX  X        X       XXX  XXX  XX       X  XXXX  XX       X    XX 

X  XX   X   X  XXX  XXX  XXXXX       XX 
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k  .   RECONNECTION 

h .1  General 

Another  form  of  noise  is  burst  noise.   It  differs  from  random 
noise  because  it  is  completely  local;  all  points  in  a  small  area  have  been 
erased.   Clearly,  the  noise  cleaning  routine  cannot  be  used  to  correct  this 
noise.   The  method  of  correction  is  to  reconnect  the  lines  which  have  been 
broken. 

To  select  the  lines  which  are  to  be  reconnected,  one  must  consider 
such  things  as  the  direction  and  the  curvature  of  the  lines,  and  the  size  of 
the  gaps.   It  is  important  that  the  decision  consider  the  continuity  of  the 
curvature  of  the  lines . 

k  .2     Generation  of  Burst  Noise 

The  pattern  used  in  the  noise  cleaning  testing  is  also  used  in  the 
testing  of  the  reconnection  routine.   The  point  where  the  curved  line  in  Figure  3 
crosses  itself  was  selected  for  the  center  of  the  noise  area. 

A  plane  called  NOISE  is  used  as  a  negative  mask  on  the  input  picture. 
Each  time  the  program  is  iterated,  the  points  in  NOISE  are  smeared  to  produce 
a  larger  noise  area.   Then  a  plane  called  NOISY  is  set  as  NOISY  =  ORIG  •  NOISE, 
thus  causing  points  to  be  removed  from  ORIG. 

k  .3  Reconnection  Process 

The  first  step  in  the  reconnection  process  is  the  selection  of  end 
points.   This  is  done  by  nearly  the  same  process  used  in  the  noise  cleaning 
routine.   The  difference  is  that  the  labels  are  extended  to  nonlabeled  points 
after  the  LABEL  routine.   The  nonlabled  points  are  given  the  labels  of  their 
neighbors  by  a  directional  smearing  process. 

The  end  point  coordinates  are  listed,  then  the  program  leaves  the 
parallel  processing  mode  and  processes  this  list. 


1.   See  program  listing  in  Appendix  II  and  flow  chart  in  Figure  k, 
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Read  ORIG  picture 


i 


Add  NOISE 


1 


LABEL  all  points 'of  NOISY  picture 


I 


Find  END  POINTS  in  picture 


1 


List  and  Count  END  POINTS 


I 


For  each  pair  of  END  POINTS: 
Test  for  possible  connections  and 
calculate  continuity  of  curvature  through  gap 
(see  Figure  ^b) 


I 


J 


For  each  END  POINT; 
Find  the  best  connection  point  (if  one  exists) 
by  selecting  minimum  change  in  curvature  through  the  gap 


I 


Reconnect  the  best  connections 


I 


Output:   Reconnected  picture 


Figure  ka .      Reconnection  Routine  Flow  Chart 
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Enter  from  Fig.  Ua 


Get  next  pair  of  END  POINTS 


1 


Return  to  Fig.  ka 


(No  more  pairs) 


L^  Length  of  gap  between  END  POINTS 


Record  BAD-PAIR 
Indicator 


Find  three  segments  of  length  L  in  each  line 


1 


Find  the  directions  and  weighted 
dire ctions  of  the  s e gme n t s 


I 


Find  curvatures  of  gap  and  outside  sections 
(using  unweighted  directions) 


Does  the  sign  of  the  gap  curvature  agree 
with  those  of  the  outside  sections? 


I 


No 


Yes 


Find  curvatures  on  both  sides  of  gap 
(using  weighted  directions) 


I 


Put  differences  between  curvatures  in  ' DIFERS'  array 


Figure  l*b .   END  POINT  Test  Flow  Chart 
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For  each  pair  cf  end  points ,  the  length  of  the  gap  which  may 
separate  them  is  approximated  as  the  larger  of  the  x  or  y  differences  in 
coordinates.   If  this  length,  L,  is  greater  than  12,  the  gap  is  too  large  for 
the  size  of  the  burst  being  used,  and  the  next  set  of  points  is  started. 

If  L  <  12,  the  subroutine  SPREAD  is  used  to  find  line  segments 
of  length  L  which  extend  from  the  end  points.   This  is  done  three  times  for 
each  end  point.   The  second  segment  begins  at  the  far  end  of  the  first  segment, 
and  the  third  at  the  end  of  the  second.  This  produces  three  sequential  segments 
of  each  line  (seven  in  all,  including  the  gap)  of  length  L0   The  eight  coordinate 
sets  of  the  junction  and  terminal  points  are  put  in  a  block  called  SEGSC  (segment 
coordinates) „ 

The  SPREAD  routine  is  entered  with  the  starting  coordinates  in 
location  ENDC  and  L  in  PAX  index  register  8„   SPREAD  returns  to  the  main  program 
with  the  other  end  point  of  the  segment  m  ENDC, 

Then,  for  each  segment.,  the  direction  of  the  line  is  found  by 
computing  DIR  =  ARCTANGENT  (DY/DX)„   This  direction  is  changed  to  conform  to 
our  labeling  convention  by  making  1,0  <  DIR  <  9°0° 

With  seven  segments,  there  are  six  junction  points  (three  on  each 
side  of  the  gap),  At  each  junction  point  the  change  in  direction,  DELTAD,  is 
computed  as  the  difference  between  the  two  segment  directions.  This  gives  an 
indication  of  the  curvature  of  the  line  at  that  point. 

Also,  a  weighted  set  of  changes  in  direction  are  computed  by 
multiplying  the  DELTAD1 s  nearest  the  center  of  the  string  by  3i    multiplying 
the  next  nearest  DELTAD' s  by  1,5,  and  by  leaving  the  outside  DELTAD' s  un- 
modified,  'These  weighted  values  are  put  in  a  block  called  WATDD,   This  gives 
more  emphasis  to  the  curvatures  near  the  gap,  and  is  useful  in  getting  an 
accurate  indication  of  the  continuity  of  the  curvature, 

A  preliminary  test  of  the  signs  of  the  curvatures  , is  made.   Three 
curvatures  are  computed:   one  is  the  sum  of  the  two  center  DELTAD' s,  the 
other  two  are  the  sums  of  the  extreme  DELTAD:s,   If  the  signs  of  these  two 
outside  curvatures  agree,  the  center  (or  gap)  curvature  is  tested „   If  its 
sign  disagrees  with  the  other  two,  computation  is  terminated  and  the  next 
pair  of  end  points  is  started.   If  the  signs  of  the  two  outside  curvatures 
disagree,  no  conclusions  can  be  drawn,  and  the  computation  continues. 
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Two  more  curvatures  are  computed.   They  are  the  sums  of  the  first 
three  and  last  three  WATDD's.   The  difference  between  these  two  is  put  in  a 
two-dimensional  array  called  DIFERS.   DIFERS(a,b)  and  DIFERS(b,a)  are  set  to 
this  difference  for  the  end  points  a  and  b.   If  the  computation  terminates 
prematurely  (L  >  12,  sign  test  fails,  etc.)  these  DIFERS  are  set  to  10.0 
(larger  than  any  possible  difference). 

This  completes  the  computations  for  each  pair  of  end  points.   When 
all  pairs  have  been  computed,  the  test  for  the  best  connections  begins. 

First,  the  value  10.0  is  stored  in  all  primary  diagonal  elements  of 
DIFERS.   Then,  for  each  value  of  a,  the  minimum  value  of  DIFERS(a,b)  is  found 
and  this  b  is  the  best  connection  for  a.   The  end  points  a  and  b  are  entered 
in  the  reconnection  list,  OKCS  (O.K.  connections).   If  no  best  connection  is 
found,  no  points  are  entered  in  OKCS. 

When  all  a's  have  been  checked,  each  pair  of  points  in  OKCS  is 
connected  in  NOISY  and  this  becomes  the  output  picture. 


h.k     Results 

The  results  of  the  program  are  shown  in  Figure  5„   Three  sizes 

of  burst  were  tried,  and  the  outputs  show  that  the  errors  in  reconnection 
increase  as  the  size  of  the  burst  increases. 

However,  these  results  show  that  the  routine  is  sound  in  principle 

and  can  be  used  in  the  cleaning  of  burst  noise  of  small  size. 
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EXTRACTION  OF  FEATURES 

■ 

.ion  01  minary  approach  to  the  extraction  of 

atures  from  a  Line-like  drawing.   The  program  has  not  been  written,  but 
would  rely  on  the  reconnection  process  of  Part  k.      The  purpose  of  the  program 
will  be  only  a  low-level  extraction,  with  no  identification  of  the  extracted 
portions.   (Further  processing  would  require  the  identification,) 


Extraction  Process 

As  an  example  of  this  extraction  process,  consider  the  face  shown 
below: 


For  the  extraction,  it  is  desirable  to  separate  the  ears,  eyes,  nose,  mouth 
and  neck  from  the  outline  of  the  face. 

The  eyes,  nose  and  mouth  are  disjointed  from  the  other  parts  and 
are  easily  extracted.   The  neck  and  ears  are  more  difficult  because  they 
are  directly  connected  to  the  face  outline. 

The  first  step  in  extracting  the  features  is  to  label  the  points 
of  discontinuity  ("punctuation"  points)  as  follows: 

node  -  lines  enter  the  point  from  four  directions;  i.e.*^C 

branch  -  lines  enter  the  point  from  three  directions;  i.e.   / — 

vertex  -  lines  enter  the  point  from  two  nonopposite  directions,  i,e.^ 

terminal  -  only  one  line  enters  the  point;  i.e.  *""" 

-31- 


Directional  continuity  through  punctuation  points  is  then  tested. 

This  results  in  the  extraction  of  the  neck  and  the  ears,  but  leaves  the  outline 

of  the  face  in  two  segments.   The  reconnection  routine  is  then  used  to  correct 
such  errors o 

In  more  difficult  cases,  higher  level  tests  will  be  necessary. 
These  could  include  a  combination  of  the  direction  continuity  test  and  the 
reconnection  probability,  or  could  make  use  of  a  dictionary  of  basic  forms 
found  in  the  picture  [l]. 


[1]   Ledley,  R.  S.  and  J,  B,  Wilson,  "Concept  Analysis  by  Syntax  Processing." 

Proc,  of  Amer.  Doc.  Inst.,  Annual  Meeting,  Vol.  1,  pp.  1-8  (October,  196U) 
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APPENDIX  I.   PLOTXY  ROUTINE 

An  addition  to  PAX  has  been  used  in  these  programs.   The  addition 
is  called  PLOTXY.   Given  a  series  of  coordinates,  this  routine  plots  through 
the  series  by  connecting  the  coordinates  with  straight  line  segments. 

The  routine  is  called  by  the  order: 

PLOTXY     PLANE, PX,FIRSTC,CX,N, NX 

where  PLANE, PX  is  the  plane  (indexed)  where  the  line  is  to  be  plotted, 

FIRST, CX  is  the  location  of  the  first  coordinate  (indexed)  to  be  plotted, 
and   N,NX     is  the  number  of  coordinate  sets  (indexed)  in  the  list. 

The  previous  contents  of  the  plane  are  undisturbed. 
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APPENDIX  II      ♦**      PROGRAM  LISTINGS 
LABEL  SUBROUTINE      -   PAGE  1  OF  1   - 


LABEL 

INDEX 

SET,  1,4 

BCOLOP 

NEWVAL, , CLEAR 

LOOPA 

INDEX 

DECK,  1, 1 

BOOLOP 

TEMPI, , EQUAL, INPUT 

BCOLOP 

TEMP2, , EQUAL, INPUT 

SHIFTD 

TEMPI,, 1, 1,2 

SHIFTD 

TEMP2, ,4, 1,2 

BCOLOP 

TEMPI, , AND, TEMPI, ,TEMP2 

BOOLOP 

IODIR, 1, AND, TEMPI, , INPUT 

BCOLOP 

I0DIR+4*144, 1, EQUAL, I0D1R,1 

BCOLOP 

NEWVAL, , OR, NEWVAL, , IODIR, 1 

JUMP 

NOZERO, LOOPA, 1 

INDEX 

SET, 1,4 

BOOLOP 

NOCON, , EQUAL, NEW VAL,, , ,BAR 

LOOPU 

INDEX 

DECR, 1,  1 

INDEX 

TRANS, 2, 1 

INDEX 

TRANS, 3, 1 

INDEX 

TRANS, 4, 1 

INDEX 

INCk,2,4 

INDEX 

INCR,3, 3 

INDEX 

0ECK,4, 1 

TMARK 

DL1,  1, TEMPI, , NOCON, ,  IODIR, 2, EQ,1 

INDEX 

SET,  11,3 

LOOPU1 

INDEX 

DECR, 11,  1 

JUMP 

LESS, GOON1A, 3,8 

INDEX 

DECR, 3, 8 

GOON1A 

JUMP 

MORE, GOON1B, 4,-1 

INOEX 

INCR,4,8 

GOON1B 

JUMP 

LESS, GOON1C, 4,8 

INDEX 

DECR, 4, 8 

GOON1C 

TMARK 

DL1,3,TEMP2, , TEMPI,, IODIR, 4, EQ,1 

GOON2 

BOOLOP 

NEWVAL, , OR, NEWVAL, , TEMP 2 

BCOLOP 

NEWVAL, , OR, NEWVAL, , TEMP2 

BCOLOP 

IODIR, 1, OR, IODIR, 1, TEMP 2 

BCOLOP 

IODIR, 3, OR, IODIR, 3, TEMP 2 

INDEX 

INCR,3, 1 

INDEX 

INCR,4,1 

JUMP 

NOZERO, L00P81, 11 

JUMP 

NOZERO, LOOPB,  1 

BOOLOP 

NOCON, , EQUAL, NEW VAL,, ,,BAR 

JUMP 

RETURN, ,8 

CLA 

DL 

2,1,5 

DL 

2,2,6 

DL 

2,3,7 

DL 

2,4,8 
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NOISE  CLEANING  ROUTINE      -   PAGE  1  OF  Z       - 

SETMOD  11 

BCOLOP  ALL,  , CLEAR, ,,, ,BAR 

NEWPIC  r<E4D  ORIG 

LETTER  (  X) 

PRINT  ORIG, ,72, 72 

COM^EN  (-ORIGINAL  PICTURE) 

TMARK  CL4LL, , INPUT, , ALL, , OR IG, ,GE , 5      LABEL  T3 

JUMP  LINK, LABEL, 8 

INDEX  SET, 1,4 

BCOLOP  T5ALL, , EQUAL ,NEWVAL 

600LOP  T5BAR, , EQUAL ,NOCON 

LOOPl   INDEX  DECK,  1,1 

BCOLOP  T5, 1 ,OR,  IOO IK, 1,  IOD IR+4» 144, 1 

JUMP  N0Z6R0, LOOPl, 1 

TMARK  GLAlL,  ,  INPUT  ,  ,  ALL  ,  ,  ORIG,  ,  GE  ,4       LABEL  T4 

JUMP  LIN.<,LABEL,8 

INDEX  SET,  1,4 

L00P2   INDEX  DECR,1,1 

BOO LOP  T4,  1,0R, IODIR, 1,  I OD IK +  4* 144  ,  1 

JUMP  NOZERO, LOOP?, 1 

EXTEND  INDEX  SET, 1,3        BEGIN  EXTENSION  OF  T5  INTO  T4 

INDEX  SET, 2,4 

LABELE  INDEX  DECR,1,1       LABEL  END  POINTS 

INDEX  DECK, 2,1 

JUMP  MORE, G00N3, 2,-1 

INDEX  INCR,2,4 

GCCN3   TMARK  DLC , 1 , TEMPI , , T5, 2 , T 5B AK , , LQ , 3 

BCOLOP  TEM<>2f  »  EQUAL,  T5ALL 

SHIFTD  TEMP2,, 1,1,1 

BCOLOP  CENDS, 1, AND, TEMPI, , TEMP2 

JUMP  NOZERO, LABELE, 1 

INDEX  SET, 1,3 

XLA3S   INOEX  DECK, 1,1       EXTEND  LAbELS  RY  ADDING  T4 ( D- 1 ) + T4 ( 0+ 1 ) 

INDEX  TRANS, 2,1 

INDEX  TRAslS,3,l 

INDEX  DECK, 2,1 

INDEX  INCK,3,1 

JUMP  MORE, G00N4A, 2,-1 

INDEX  INCR,2,4 

GOOiN-tA  JUMP  M0Rc,G00N4B,  3,-1 

INDEX  INCR,3,4 

G0CN4B  JUMP  LESS, G00N4C, 2,4 

INDEX  DECK, 2, 4 

G0CN4C  JUMP  LESS, G00N4D, 3,4 

INDEX  DECR,3,4 

GC0N4C  BCOLOP  DENOSX , 1 , OR , DENDS , 1 , T4 , 2 

BOOLOP  DENDSX, 1, OR, DENOSX, 1, T4, 3 

JUMP  NOZERO, XLABS, 1 
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NOISE  CLEANING  ROUTINE      -   PAGE  2  OF  2   - 

INDEX    SET, 1,8 

INDEX    SET, 2, 4 
GROW    INDEX    DECR,1,1       EXTEND  DENDSX  INTO  T4 

INDEX    DECR,2,1 

JUMP     MORE, GU0N6, 2,-1 

INDEX    INCR,2,4 
GCCNo   CHAIN    DL1,  1,NEWPTS,1, DENDSX, 1, T4, 2 

JUMP     NOZERO,GROW,  1 

INDEX    SET, 1,3 

INDEX   SET, 2, 4 
COMP'XR  INDEX    DECK, 1,1       WERE  NEW  POINTS  PICKED  UP 

INDEX    DECR,2,1 

JUMP     MORE, GOON5, 2,-1 

INDEX    INCR,2,4 
GC0N5   BCOLOP   COMP , , AND , NE wP  T S ,  1 , T5 , 2 , BAR 

JUMP     NONULL,'<E  ITER,  COMP 

JUMP     NOZ£RO,COMPAR, i 

LETTER   (  E A1N23+345-6*. %)  NO,  PRINT  RESULTS 

CPRINT   4,72, 72, T5,, T 5+144, ,T 5+2* 144, , T 5+3*1 44 

COMMEN   (-PROCESSED  PICTURE) 

JUMP     UCNEWPIC 
REITr.R  INDEX    SET,1,4  YES,  REITERATE  EXTENSION 

ADOPTS  INDEX    DECRiltl 

BCOLOP   TEMPI,  ,OR,NEWPTS, 1 , NE WPTS+4* 144 ,  1 

BOOLOP   T5, l,OR,l 5, 1, TEMPI 

BCOLOP   T5BAR, , AND, T5BAR, .TEMPI, , BAR 

JUMP     NO/.ERO,  ADOPTS,  1 

JUMP     UC, EXTEND 
CLC     DL     3,8,1,2 

DL     3,1,2,3 

OL     3,2,3,4 

DL     3,3,4,5 

DL     3,4,5,6 

DL     3,5,6,7 

DL     3,6,7,8 

DL     3,7,8,1 
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REC3NNECTI0N  ROUTINE      -   PAG  =  1  OF  3 

SETMOD  72 

READ  OR  I  3 

BCULOP  NOiSEt tCLEAR 

WRITE  NOISE, ,21, ,48 

INDEX  SET,  12,1 

ADDNCS  INDEX  INCH, 12,1  ADD  NOISE 

JUMP  MORE, DONE,  12»  5  NOT  100  MUCH  NOISE 

MARK  DLALL, » TEMPI, , NOISE 

BCOLOP  NOISE, , EQUAL, TEMPI 

RCOlOP  NQISEYt *  AND, ORIC, , NOISE, , BAR 

INPUT   EQU  NOISEY  US:  NOIStY  FOR  INPU1 

JUMP  LINK, LABEL, 8  LABEL  I  ^PUR  PICTURE 

TEST1   BCOLOP  EXTRAS ,, AND , NOISEY, , NOCON     FIND  UNLABELED  POINTS 

JUMP  NULL, G00N3, EXTRAS  IF  NUNE  GO  ON  I U  3 

INDEX  SET, 1,3 

EXTL'^B  INDEX  D£CR,1,1  EXTEND  LABELS  FROM 

MARK  CLC, 1, TEMPI, , IODIR, 1, EXTRAS      NEARBY  POINTS 

BCOLOP  IODIR, 1, OR, IODIR, 1, FEMPi 

BCOLOP  NOCON, , AND, NOCON, , TEH PI,, BAR 

JUMP  N0ZuRO,cXl'LAB,  1 

JUMP  UCTEST1 

GCONJ   INDE<  SET, 1,8 

BCOLOP  ENDS,, CLEAR 

FINDES  INDEX  DECK, 1,1  FIND  END  POINTS 

TMARK  DLC1,  TEMPI,,  IODIR,  1,  NOCON,  ,  CO,  3 

BCOLOP  TEMP2, , EQUAL, IODIR,  1 

SHIFTD  TEMP2,, 1,1,1 

BCOLOP  ENDPTS,l,AND,rEMPl,,TEMP2 

BCOLOP  ENDS, , OR, ENDS, ,ENDPTS,  1 

JUMP  MOZERO, FINDES, 1 

INDEX  SET, 1,3 

LISTXY  ENDS, ,ELIST,  15  LIST  END  POINTS 

JUMP  M0RE,G00F1,21, 15  TEST  FOR  TOO  MANY 

JUMP  LESS, G00F4, 21,2  TOO  FEW 

PAX  0,7  SET  UP  PARAMETERS  IN  PROGRAM 

SXA  TESTC-2,7 

TXI  •♦1,7,-1 

SXD  G00N13+1,7 

ALS  18 

STD  C0UNT1+1 

SUB  DECR2 

STD  C0UNT2+1 

AXC  »  +  2,7 

TRA  PAX 

INDEX  SET, 6,0 

INDEX  SEF,7,0 

NXTPT   INDEX  IlMCR,7,l  DO  FOR  EVERY  PAIR 

C0UNT1  JUMP  LESS,G00N4,7,»*      *»  =  NO.  OF  ENDS 

INDEX  INCR.6,1 

C0UNr2  JUMP  MOKE, A0DLIN,6,*»         **  =  NO.  OF  ENDS  -2 

INDEX  TRANS, 7, 6 

INDEX  INC.K,7,1 
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RECONNECT  ION  ROUUNt 


pa;;:  2  of  a 


GCCN-t 


J7094C 

LAC 

LAC 

CLA 

STO 

ARS 

STO 

CLA 

STO 

ARS 

SUB 

STO 

SLrt 

STZ 

STZ 

CLA 

STA 

CLA 

STA 

CLA 

SUB 

STO 

SSP 

CAS 

TRA 

TRA 

CLA 

STO 

CAS 

TKA 

NOP 

AXC 

TRA 

BCOLOP 

J7094C 

CLA 

STO 

AXC 

TRA 

JUMP 

J7094C 

CL4 

STO 

AXC 

TRA 

JUMP 

J7094C 

CLA 

STO 

AXC 


*  +  l 

IR- 1+6,1 

IR-1+7,2 

ELIST,  1 

SEGCS+3 

13 

TEMPI 

EL  I  ST, 2 

SEGCS+-4 

18 

TEMPI 

DX 

MAGOX 

TEMPI 

TEMP2 

ELIST, 2 

TEMPI 

EL  I  ST,  1 

TEMP2 

TEMP2 

TEMPI 

DY 

MAGDX 

•  ♦3 
»+2 

MAGDX 
IR-1+8 
FX12 
G00F2 


GET  THE  PAIR  OF  COORDINATES 
FIND  LENGTH  =  MAX (MAGDX ,MAGOY) 


IR8  =  L 


AL, ,EQUAL,NOCON,, ,,BAR 

SPREAD  SEGMENTS  OF  LENGTH  L 
S*3  INTO  LINES  FROM  ENDS 


*  +  2, 

PAX 

NEW'/ 

•♦1 

SEGC 

ENDC 

*  +  2,7 
PAX 
LINK,SPR6A0,9 

*  +  l 
ENDC 
SEG 

*  +  2, 
RAX 

LINK, SPREAD, 9 

»  +  l 

ENDC 

SEGC 

*  +  2, 


S  +  2 

7 


S  +  l 

7 
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FINCD 


TRA 

PAX 

JUMP 

LINK, SPREAD, 9 

J7094C 

**1 

CLA 

ENDC 

STO 

SEGCS 

CLA 

SEGCS+4 

STO 

ENDC 

AXC 

*  +  2,7 

TRA 

PAX 

JUMP 

LINK, SPREAD, 9 

J7094C 

•  +  1 

CLA 

ENDC 

STO 

SEGCS+5 

AXC 

*  +  2,7 

TRA 

PAX 

JUMP 

LlNK,SP« 

EAO,  9 

J7094C 

»  +  l 

CLA 

ENDC 

STO 

SEGCS+6 

AXC 

*  +  2,7 

TRA 

PAX 

JUMP 

LINK, SPREAD, 9 

J7094C 

»  +  i 

CLA 

ENDC 

STO 

SEGCS+7 

AXT 

7,1 

CLA 

SEGCS+7, 

1 

ARS 

19 

STO 

TEMPI 

CLA 

SEGCS+8, 

I 

ARS 

18 

SUB 

TEMPI 

STO 

DX 

STZ 

TEMP2 

STZ 

TEMPI 

CLA 

SEGCS+7, 

I 

STA 

TEMPI 

CLA 

SEGCS+8, 

1 

STA 

TEMP2 

CLA 

TEMP2 

SUB 

TEMPI 

STO 

DY 

CLA 

FXTOFL 

ORS 

DX 

ORS 

DY 

ZAC 

FAO 

DX 

STO 

DX 

ZAC 

FAD 

DY 

FIND  DIRECT IOM  OF  SEGMENTS 


CONVERT  DX  AND  DY  TO  FLOATING  POINT 
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RECONNECT  ION  ROUTINE 


PAG-  4  OF  8 


STO 

OY 

NZT 

OX 

TRA 

NOOX 

Nzr 

OY 

TRA 

NOD/ 

CLA 

OY 

FDP 

OX 

STQ 

OIR*7,  1 

TSX 

ATAN,4 

TXH 

DIR+7,  1 

STO 

D I R  +  7 ,  1 

CLA 

DX 

TMI 

NEGDX 

CLA 

DY 

TMI 

NEGOY 

LDO 

DIR  +  7,  1 

G00N9 

FMP 

F40VPI 

FAD 

FL1 

CAS 

FL9 

NOP 

FSB 

FL8 

CAS 

FL1 

TRA 

♦  +  3 

TRA 

*  +  2 

FAD 

FL8 

GC0N8 

STO 

D I R  +  7  , 1 

TIX 

FINDD,1, 1 

AXT 

6,1 

CURVG 

CLA 

0IR+6, 1 

FSB 

DIR+7, 1 

STO 

OELTAD+6,  1 

SSP 

CAS 

FL4 

TRA 

HIGHDD 

NOP 

GCCNIO 

CLA 

XCA 

DELf AD  +  6,  1 

FMP 

WATES+6, i 

STO 

WATOD  +  6,  I 

TIX 

CURVS, 1,  1 

STEST 

CLA 

DELTAD 

FAD 

DELTAD+1 

STO 

COMMON 

CLA 

DELIAD+2 

FAD 

DELTAD+3 

STO 

COMMON+1 

CLA 

DELTAD+4 

FAD 

DELI AD+5 

TPL 

PTEST 

CLA 

COMMON 

AT  AN  =  SYSTEM  ARCTANGENT  ROUTINE 

EXI I S  WITH  -PI/2  .Lc.  IHE1 A 
DIRECT  ION  =  ATAN( DY/OX) 
JIUL1    CONVERT  TO  l.LE.  0 1 R  .LT.'i 


.Li 


PI/? 


FIND  CURVATURES  AMD  WEIGHTED  CURVATU'.ltS 


TEST  CONTINUITY  OF 

SIGN  OF  CURVATURE 
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RECJNNECTICN  ROUTINE 


PAG..  5  OF  8 


GCCNll 


PTEST 


NODX 


D7 


NODY 


C5 


NEGDX 


NEGDY 


HIGHCD 


TPL 

GOON  11 

CLA 

COMMON +1 

FPL 

GOOF2 

ZAC 

AXT 

3tl 

PAD 

wATDD+3,  1 

nx 

♦-1,1,1 

STO 

CDAV1 

ZAC 

AXT 

3,1 

FAD 

WATOD+6,  1 

TIX 

*-l,  1,1 

FSB 

cda;i 

SLW 

CCD 

TSX 

SETXS,4 

CLA 

COD 

STO 

OIFtRS, 7 

STO 

DIFERS.6 

AKC 

*  +  2,7 

TRA 

PAX 

JUMP 

UCtNXTPT 

CLA 

COMMON 

TMI 

GOON11 

CLA 

COMMON+1 

TPL 

GOON  11 

TRA 

G0OF2 

CLA 

DY 

TMI 

D7 

CLA 

FL3 

TRA 

G0ON8 

CLA 

FL7 

TRA 

GCON8 

CLA 

DX 

TMI 

05 

CLA 

FL1 

TRA 

GC0N8 

CLA 

FL5 

TRA 

GC0N8 

CLA 

OIRi-7,1 

FAD 

PI 

XCA 

TRA 

GOON  9 

CLA 

DIR+7,1 

FAD 

TWOPI 

XCA 

TRA 

GOON  9 

CLA 

DELTAD+6, 1 

TMI 

*  +  3 

FSB 

FL8 

TRA 

•  ♦2 

FIND  AVERAGE  EFFECIIVE  CUKVATU^c  AS 
SUM  OF  WEIGHTED  CHANGES 
IN  C  IKECT  ION 


FIND  DIFFERENCE  IN  EFFECTIVE 
CURVATURES  AND  PUT  IN 
'DIFERS'  ARRAY  F0.<  FUTURE 


TEST 
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FAD  FL8 

STO  DELTAD+6,1 

TRA  GOOMIO 
ADOLIN  J7094C   *+l  PREPARE  TO  ADD  LINES  TO  ORIG 

AXT  0,1 

CLA  FLU  STORE  ERRUR  INDICATOR  (10.0)  ON  PRIMARY 

AXT  15,2  DIAGONAL  OF  DIFERS  ARRAY 

STO  DIFERS, 1 

TXI  *+l,l,l6 

TIX  «-2,2,l 

AXT  0,5 

AXT  0,1 

NXTPNT  PXA  0,1 

XCA 

MPY  FX15 

XCA 

PAX  0,2 

CLA  FL9P6 

TESTC   CLA  DIFERS, 2       TEST  FOR  BEST  CONNECTION 

CAS  MIND 

TRA  GC0N12 

NOP 

STO  MIND 

SXA  BESTC4 

G0CN12  TXI  *+l,4,l 

TXI  *+l,2,l 

TIX  TESTC, 3,1 

CLA  MIND 

CAS  FL9P3 

TRA  G00N13       NO  POSSIBILITIES 

NOP 

LAC  6ESTC6  ENTER     IN    LIST    OF    O.K.     COORDINATES 

CLA  ELIST,6 

STO  0KCS,5 

TXI  »+l,5,-l 

SXA  BESTCfl 

LAC  BESTC6 

CLA  ELIST.6 

STO  CKCS,5 

TXI  *+l,5,-l 

G0GN13  TXI  ft+lfltl 

TXL  NXTPNT, i,*«    **  =  ENDPTS  -  1 

PCA  0,5 

TZE  NOCONN         IF  NO  CONNECTIONS  ARE  FOUND 

ARS  1 

STO  IR-l+2 

AXC  *+2,7 

TRA  PAX 

INDEX  SET, 1,0 
ADOPTS  PLOTXY   NOI SEY , , OKCS , 1 , 2       ADD  POINTS  TO  RECONNECT 
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RECONNECTION  ROUTINE 
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SETXS 


CONE 
SPREAD 


GROW 


TESTES 


GOCN7 
GOOF1 

G00F3 
G0QF2 


GCCF4 
NCCCNN 


INDEX 

INDEX 

JUMP 

LETTER 

PRINT 

JUMP 

LDQ 

MPY 

XCA 

ADD 

PAX 

LDQ 

MPY 

XCA 

ACD 

PAX 

TRA 

HALT 

INDEX 

BCOLOP 

WRITEP 

BCOLOP 

INDEX 

MARK 

BOOLOP 

BCOLOP 

JUMP 

LISTXY 

JUMP 

JUMP 

BCOFUN 

JUMP 

JUMP 

COMMEN 

JUMP 

J7094C 

TSX 

CLA 

STO 

STO 

AXC 

TRA 

JUMP 

COMMEN 

JUMP 

AXC 

TRA 

COMMEN 

JUMP 


DECR,2,  1 

INCR,  1,2 

MORE, ADOPTS, 2,0 

(  X) 

NOISEY, ,72,72 

UC, ADDNOS 

IR-1+6  SETS 

FX15 

IR-1+7 
0,7 

IR-1+7 
FX15 


7094  INDEX  REGISTERS  TO  INDEX 
2  -  DIMENSIONAL  DIFERS  ARRAY 


IR- 

0,6 

1,4 

DUM 

TRA 

END 

END 

NEW 

DEC 

DLS 

NEW 

END 

NOZ 

END 

EXA 

LES 

END 

UC, 

RET 

(-. 

UC, 

•♦1 

SET 

FL1 

DIF 

DIF 

♦  ♦2 

PAX 

UC, 

(- 

UC, 

«  +  2 

PAX 

{- 

UC, 


1+6 


p  **   AT  LAST   ** 

NS,  1, 8 

, .CLEAR 

, ,ENDC 

VAL, , AND,NEWVAL, , END, , BAR 

R,i,l       GROW  L  POINTS  BACK 

UR, ,1 EMP1, ,END, ,NEWVAL 

VAL, ,  AND,NEWVAL, , TEMPI , ,BAR 

, , EQUAL, TEMPI 

ERO,GROW, 1 

, ,ENDC,2 


ALONG  LINE 


CT,G00F3, 21,0    EXIT 
S,GOON7,21,3     TEST 


IF  WOULDN'T  GROW 
FOR  MULTIPLE  ENDS 

REMOVE 


, ,END, , END,, ( 13+15+17+35+37+57$) 

TESTES 

URN, ,9 

GT.  15 

ADDNOS 


SOME 


ENC  POINTS) 


XS,4 
0 

ERS,6 
ERS,  7 
,7 

NXTPT 
.LT.  2 

ADDNOS 
,7 


PUT  'BAD  SET  OF  ENDS' 
IN  DIFERS  ARRAY 


INDICATOR 


END  POINTS) 


NO  CONNECTIONS) 
ADDNOS 
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CLC 


DECR2 

FXTOFL 

FL1 

FL3 

FLA 

FL5 

FL7 

FL8 

FL9 

FL9P3 

FL9P6 

FLIO 

FX12 

FX15 

WATES 

DIFEKS 

SEGCS 

CKCS 

COMMCN 

ELIST 

ENDC 

DIR 

DELTAD 

WATDD 


RECONNECTION  ROUTINE 

CONSTANTS  AND  BLOCK  RESERVATIONS 

DL  3,8,1,2 

DL  3,1,2,3 

OL  3,2,3,4 

OL  3,3,4,5 

DL  3,4,5,6 

DL  3,5,6,7 

DL  3,6,7,8 

DL  3,7,8,1 

OCT      2000000 

OCT      233000000000 

DEC      1.0 

DEC      3.0 

DEC      4.0 

DEC      5.0 

DEC      7.0 

DEC      8.0 

DEC      9.0 

DEC      9.3 

DEC      9.6 

DEC      10.0 

DEC      12 

DEC      15 

ASSIGN   DX,0Y,MAGDX,MIND,BESTC,DDAV1,DD0 

DEC      1.0,1.5,3.0,3.0,1.5,1.0 

BES      15«15 
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BES 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 
BSS 


15*L5 

8 

30 

3 

16 

3 

7 

6 

6 
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